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ANALYSIS OF DEVELOPING LAMINAE FLOW AND HEAT TRANSFER IN A TUBE FOR A GAS WITH VARIABLE PROPERTIES 


by Robert G. Deissler and Alden F. Presler 
Levis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

Laminar tube flow and heat transfer for helium gas are analyzed numerically by using the compressible 
Navier-Stokes, energy, and continuity equations in finite difference form. Both radial and axial property 
variations are considered. The radial as veil as the axial velocity are retained in the equations; the 
only assumptions made in the equations are the usual boundary- layer assumptions. The heat flux at the vail 
and the initial velocity and temperature profiles are taken to be uniform. The effects of both heat flux 
and Mach number on the Nusselt number and friction correlations are obtained. 


AUSZUG 

Laminar stromung in einer Rohre und Warmeubertragung fur Heliumgas wurden rait Hilfe der kompressiblen 
Forael von Navier-Stokes und der Energie- und Kontinuitatsformeln in endlicher Differenzform numerisch 
analysiert. Abweichungen der radialen und axialen Eigenschaften werden in Betracht gezogen. Sowohl die 
radiale vie die axiale Geschvindigkeit sind in die Gleichungen eingesetzt und nur die ublichen Grenzschicht- 
annahmen sind in den Gleichungen gemacht. Der Warmestrom an der Wand und die Ausgangsgeschwindigkeit und 
-temperaturprofile sind als glei chinas sig angenommen. Der Einfluss des Warmestroms und der Machzahl auf die 
Nusseltzahl und die Reibungsbeziehungen werden erhalten. 


AHHOTAUHfl 

MucjieHHbiM cnocoGoM ajih reran zccjreflyioTCii jiaMHHapHUH noTOK b Tpyfixe h Teruionepeflaua 
npHMeHHH ypaBHeHZH HaBbe-CTOKca, OHeprHH n Hepa3puBH0CTn ajih c*HMaeMOro ra3a b BH^e KOHeu- 
hlix pa3HOCTeH. PaccMaTpHBaioTCH Kan pa^HajibHbie Tax h oceBtie nepeMeHHbie CBOHCTBa; b ypaB- 
HeHunx cflejiaHbi tojibko ofeiuHbie npefliiojiosceHH.fi OTHOCHTejibHO norpaHHunoro cjioh. IIpeflnojiaraeTCH 
paBHOMepHOCTb nOTOKa Tenjia na CTeHKe h npo$njieH HauajibHOH cicopoeTH h TeMnepaTypbi. IIojiyueHbi 
bjihhhhh nOTOKa Temia h uzcjia Maxa Ha uhcjio HyccejibTa h Ha Koppejimizu TpeHZH. 


INTRODUCTION 

The laminar flow and heat transfer in the en- 
trance regions of passages have been extensively 
analyzed by a number of workers. The constant- 
property Graetz problem [l and 2],1 in which the 
velocity profile is assumed fully developed, and 
where the temperature profile is calculated as a 
function of distance from the entrance, has been 
exhaustively explored. Whenever the development of 
the velocity profile has been considered, integral 
methods have generally been employed [3], or the 
nonlinear acceleration terms in the equations of 
motion have been approximated by linear terms 
[4 and 5]. A numerical solution for constant- 
property plane -Poiseuille flow in which the non- 
linear terms are retained is given in [6]. Numeri- 
cal solutions for heat transfer with developing 
velocity and temperature profiles in a tube are 
given in [7] and [8]. In those analyses the veloc- 
ity profiles are obtained from the linearized solu- 
tion in [4] . 

For the case of laminar flow with variable 
fluid properties, only a limited amount of work has 
been done. An analysis in [9] considers flow in a 
tube far from the entrance. Radial variations of 
fluid properties are considered, but axial varia- 
tions are neglected. The effects of radial veloc- 
ity are also neglected. In [10] various radial 
flows are arbitrarily introduced into an analysis 
in an attempt to study the effects of radial flow 
on heat transfer and fricti on with variable pro p- 

lumbers in brackets denote references. 


erties. A comparison of the results with the 
experiments in [10] and [ll] indicates that the 
radial flow could have an important effect. Addi- 
tional partial solutions, which study effects of 
variable fluid properties, are given in [12], 

[13], and [14], 

In order to eliminate, as far as possible, 
the dubious assumptions made in previous analyses, 
the equations of motion and energy for a gas with 
variable properties are herein solved numerically. 
Only the usual boundary -layer assumptions are 
made. Those assumptions have been made in nearly 
all of the analyses for flow and heat transfer 
in tubes, and they generally give good results, 
except in the region very close to the entrance. 
The pertinent equations and their solutions for 
uniform initial velocity and temperature profiles 
in a tube with uniform wall heating will be con- 
sidered in the next section. 


ANALYSIS 

The equations of motion, energy, and state for 
a perfect gas, for a steady-state, axially- 
symmetric flow, without swirl and body forces, can 
be written as [15] 
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If it is assumed that the boundary layer thickness 
5 « z, then Equations (l), (2) and (4) reduce to 
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where the axial density gradient in Equation (3) 

(6) was eliminated by differentiating the equation of 

state (Eq. (5)). The subscript i refers to values 
at the tube inlet. The fact that v r = 0 at the 
tube axis was used in obtaining Equation (8). 

To obtain an expression for the axial pressure 
/ 7 \ gradient, which is independent of r, Equation (6) 


at the wall is evaluated: 


where, from Equation (2), p and dp/dz are essen- 
tially independent of r. The boundary -layer as- 
sumption should be valid except, possibly, for the 
region very close to the entrance. 

An expression for the radial velocity v p can 
be obtained by integrating the continuity Equa- 
tion (3) as follows: 




By using Equations (5), (6), (7), (8), and (9), 
the final set of equations to be solved can be writ- 
ten in dimensionless form as 


Av z „ dv z _ 1 - n'T- dr i f . fl - ynfp' b^ + q , Vg(l - rM?p') aj. _ ^2 v z bp']^ 
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Az 1 bz * Re.^ 



^T2 + ^ ^7 + r ^7^tJ 


( 12 ) 


In order to compare results with those in [ 9] , the 
following values are used: 

/ v 0.68 

d' = k’ = = (1 - q'T') 0 ’ 68 (13) 

c* = 1, (14a) 

P 

and 

Pr i = 0.667 (I4h) 

These properties also correspond closely to those 
for helium, so that these results should apply to 
that gas. 

To carry out the numerical solution the de- 
rivatives in Equations (10), (ll), and (12 ) are 
replaced by their finite difference forms. For 
instance, bv^/br' is replaced by 

( v z,iHAr - v i,r^r) /2Ar '> “ 1 by 

( V z,r-Ar “ ^ZjX + v z,r-*Ar^^ r ' l n the 
problem considered here, the velocity and temper- 
ature profiles are uniform at the tube inlet, and 
the heat flux at the tube wall is also uniform. 
Thus, the initial conditions can be written in di- 
mensionless form as follows: At the tube inlet 

(z* - 0), T' = 0, p* = 0, and = 1 (except at 
r' = 1, where = 0). For the boundary condi- 
tions at the wall (r' = l), v£ = 0 and 

c^T’/dr* = -l/k*. At the tube axis (r* = 0), 
chr^/dr' = 0 and bH'/br' = 0. Since the wall heat 
flux is uniform, the parameter q* is set equal to 
a constant for each numerical run. 

With the foregoing initial and boundary con- 
ditions, Equations (10), (ll), (12), (13), and (14) 
can be used in finite difference form to calculate 
v^, T', and p* at various increments of z'. 

The values of bv^/bz 1 and dT'/dz’ occurring on 
the right sides of the equations are generally not 
known at z* = 0, but they may be found by itera- 
tion by assuming values for them and then calcu- 
lating new values from the equations. For later 
Az steps, values of dv£/az ! and bx'/bz' from 
the preceding step can be substituted in the right 
sides of the equations. These might again be iter- 
ated, although it may not be necessary for suffi- 
ciently small values of Az. The calculations can 
also be carried out by a matrix method which avoids 
the use of iteration. In the calculations, the 
ratio of Az* to (Ar*) 2 must be kept sufficiently 
small to ensure stability of the solution. The 
quantity Ar' must also be made small enough so 
that cutting it in half does not change the results 
appreciably. The calculations were carried out on 
a high speed digital computing machine. With the 
dimensionless velocity and temperature profiles, 
as well as the pressure, calculated as functions 
of z 1 , R, M^, and q', the following integral 
quantities can be obtained: 


v 


z,b 



dr * 


(15) 


T b = 


T ' 
t,b 


/ l v£r’ 

T . dr ' 

1 - 

f 1 — dr' 

Jo 1 - q’T’ 

/ ! v 'r 1 

T ' 5 dr 

-t i . q-r 


r 1 

Jo 1 - q’T' 


dr 1 


where 


T ' = T ' - £ - — v * 2 

t 2 q 1 z 

r 2 -n * W 1 


Rob 


= Pb v z,b Re i (l - r^p')^^ 

^ (i . 

2 2 


0.68 


Nu, = 


k b (T ;,b - V (1 - d-T- )°- 68 (T; >t - T’ >v ) 

4(1 - 


T,b 


(1 - r^P*)v^ b Re. 


(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


where 


« _ t' 

8 b = 


, Tl) 0.68Fd] 

W 7 \^_ i / 

(22) 

\or / 

\ / w 


T 1 

(23) 

.^0.68 


v ;,b^ v; >b d . q'-^y 

f dp'/dz' . (1 - q’T^) dp'/dz' 

g,b 2(1 - TMfp')^ 


7 Z,b 


-£l 


b 2z 'p 'v 1 2 

b z,b 


(24) 


(25) 


Results calculated from these equations will be 
given in the next section. 


RESULTS AMD DISCUSSION 


Friction factor curves are presented in Fig- 
ures 1, 2, and 3, and corresponding Nusselt number 
results for heat transfer are given in Figure 4. 
Friction factors based on wall shear stress, pres- 
sure gradient, and pressure drop are shown. They 
are multiplied by Re^, since friction factor times 
Re^ is not a function of Re^ but is a function 
only of (z/D)/Reb, q’, M^, and Fr^ . The Nusselt 
number Nu^ is also a function only of those vari- 
ables . In the uniform property case, friction pa- 
rameters and Nusselt numbers approach constant 
values for large distances frcm the entrance. 

In Figures 1 to 3, the calculated friction 
factor curves for uniform properties 
(q* = Mi = 0) are compared with the linearized 
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analyses from [4] and [5], and with the experi- 
mental data from [16] . These calculated curves 
are in substantial agreement with the data and the 
linearized analyses, except for the linearized 
results for shear stress shown in Figure 1* The 
latter apparently deviate because the shear stress 
is very sensitive to the effect of the linearizing 
appr oximat i on . 

In Figure 4, the calculated Nusselt number 
curves for uniform properties are compared with the 
analysis from [8]- The main difference between the 
present analysis and that in [8] is that the latter 
is based on the linearized velocity profiles from 
[ 4] , whereas the present analysis does not use a 
linearizing approximation. The comparison shows 
that the effect of the linearizing approximation 
on the Nusselt number plot is small. Results from 

[8] for v r = 0 (not shown) are also in good agree- 
ment with those calculated here, although the orig- 
inal calculations for v r = 0 from [7] showed 
somewhat more deviation because of the larger 
increment size used there. 

Curves for moderately high heat transfer to 
the gas (q* = 13.5) and for Mj_ = 0 and = 0.06 
are also shown in Figures 1 to 4. For the runs 
shown, the maximum temperature ratio occurs 

in the entrance region and is about 2.7. At the 
point where Nu^ is a minimum, T^T^ is about 

1.5. The effect of the heat flux on the friction 
factor is much greater than that on the Nusselt 
number. Values of the friction parameters are, 
in general, increased considerably by the heat 
transfer to the gas. On the other hand, the 
Nusselt numbers in the entrance region are in- 
creased but slightly, and part of that increase 
may be due to the slightly lower Prandtl number 
used for the runs in which q’ >0. The fact that 
the effect of variable properties on the bulk 
Nusselt number correlation is small in the entrance 
region is in agreement with the experimental find- 
ing in [11]. Farther down the tube the values of 
Nusselt number dip below the fully developed value 
for uniform properties. Finally, still farther 
down the tube, the Nusselt numbers for q* > 0 
and = 0 again approach the uniform property 
value of 48/ll. This happens inasmuch as the 
absolute temperature ratio approaches 1 

because of the increase in temperature level along 
the tube, even though the heat flux is uniform. 
Thus, radially uniform fluid Tirnperties are ap- 
proached far down the tube, although the proper- 
ties still vary axially. The values of f^^Re^ 
also approach the fully developed uniform property 
value at large values of (z/D^Re-^. Therefore, the 
axial variation of properties does not seem to 
affect the values of Nu^ and of f^-^^Re^ in the 
region far frcm the entrance. The behavior of 
fg^Rek and fp^Re^ at large distances from the 
entrance is more complicated because those quanti- 
ties include pressure changes due to momentum 
effects . 

For > 0 and q* > 0 both the Nusselt 
numbers (Fig. 4 ) and shear-stress friction factors 
(Fig. 1) increase rapidly near the end of the run. 
This occurs because the local Mach number approach- 
es 1 and choking takes place. These large in- 
creases in Nusselt number and shear stress are 
evidently associated with the large accelerations 
that occur near a Mach number of 1. 

Numerical calculations of Nusselt number, 
such as those given in Figure 4, might be used in 
conjunction with measurements of heat -transfer 
coefficients in tubes to obtain the thermal con- 


ductivities of gases at high temperatures. The 
use of such calculations would allow measurements 
at reasonably high heat fluxes, so that greater 
accuracy could be obtained than would be possible 
if it were necessary to keep the heat fluxes small 
enough to approximate uniform properties. An 
advantage of this method of obtaining thermal 
conductivities at high temperatures is that radi- 
ation effects are small. 

Velocity profiles which correspond to the 
curves in Figures 1 to 4 for = 0.06 and 

q* = 13.5 are plotted in Figure 5. At z = 0 
the profile is, of course, flat by assumption. 

As z increases, the fluid near the wall accel- 
erates because of the heat addition. This causes 
a peak in the curve to develop away frcm the tube 
axis. Farther down the tube the profile tends to 
approach a fully developed parabolic shape. Then, 
as a Mach number of 1 is approached, the profile 
flattens because of the acceleration of the fluid 
near the wall. This region of flattening corre- 
sponds to the region of rising values of f Ra 
and Nujr, in Figures 1 and 4. T * 

Figures 6 and 7 ccmpare the present zero Mach 
number results for the quasi -developed region with 
those in [9], where the effect of radial variation 
of properties is considered, but axial effects and 
radial flow are neglected. Values of Nu^ and of 
the friction parameter = f^Re^ are plotted 

against temperature ratio T w /T^ for heat addition 
to the gas. The region of quasi -developed flow for 
the present results was scmewhat arbitrarily taken 
as the region beyond the minimum in the Nusselt 
number curves. 

The Nusselt number results from [9] are in 
reasonably good agreement with the present results, 
which decrease slowly with increasing values of 
T w /T^ . In the case of the friction parameter, 
although both the present results and those from 

[9] show an increase in 5^ with increasing 
T w/Tb> '& he increases shown by the present analysis 
are considerably greater. Experimental data from 

[10] also showed increases in friction considerably 
greater than those predicted in [9]. The conclu- 
sion to be drawn, then, is that axial changes and 
radial flow are not important as far as the quasi - 
developed he at -transfer correlation is concerned, 
but that those factors can have a considerable 
effect on the friction. 

SUMMARY OF RESULTS 

The calculated he at -transfer and. pres sure -dr op 
results for uniform properties agreed reasonably 
well with analyses which make use of linearizing 
approximations in the momentum equation. Linear- 
ized shear -stress results, however, showed devia- 
tion frcm the present calculations, apparently 
because the shear stress is very sensitive to the 
effect of the linearizing approximation. Calcu- 
lated pressure-gradient and pres sure -drop results 
for uniform properties were in satisfactory agree- 
ment with available experimental data. 

The variable -property friction -factor curves 
with properties evaluated at the fluid bulk temper- 
ature showed that the level of the curves can be 
increased considerably by heat transfer to the gas. 
On the other hand, the Nusselt number correlation 
based on the bulk temperature was affected but 
slightly by variable properties in the entrance 
region. Farther down the tube, the Nusselt numbers 
with large heat transfer to the gas dropped below 
the fully developed uniform property value. How- 
ever, when the Mach number approached 1 near the 
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end of the tube, both the Nusselt number and the 
friction parameter increased sharply. 

The numerical method used here for calcu- 
lating Nusselt numbers in variable -property streams 
may be useful in obtaining thermal conductivities 
of gases at high temperatures by measurement of 
heat-transfer coefficients. 

The results calculated here for the quasi - 
developed region and a Mach number of zero were 
compared with those from a previous analysis 
(NACA TN-2410) which considered the radial varia- 
tion of properties but neglected axial effects and 
radial flow. The comparison showed that axial 
effects and radial flow have but a slight effect 
on the heat transfer results, but they have a 
marked effect on the shear stress. 


SYMBOLS 

D tube diameter 

Cp specific heat at constant pressure 



c y specific heat at constant volume 
f T,b 2 ' £ w/( p b v z,b) friction factor based on wall 
shear stress with density evaluated at 
1 2 

f g,-b -■g‘[dp/d(z/D)]/(p^ ) v z ^ t) ) friction factor based 
on pressure gradient with density evalu- 
ated at Tfc 

f p,b (V 2 )(Pi “ p)/[(z-/^)pb v z,bl friction factor 
based on pressure drop with density evalu- 
ated at T-fc 

h 9/(^t,w ~ ^t b^> transfer coefficient 

k thermal conductivity of gas 

k* k/kj: 

% v Zji /V7rRTi) 

Nu^ h D/k^ Nusselt number with conductivity 

evaluated at T^ 
p pressure 

P 1 (Pj. - v)/{p i yl )i ) 

c p,iPi/ k i Prandtl number at inlet 
q heat transfer per unit area from wall to 

fluid 

r distance from tube axis 

r w tube radius 

r' r/r w 

Re i Pi v z,i ®/^i inlet Reynolds number 

R perfect gas constant 

Reb Pb v z,b C/Mb local bulk Reynolds number 
T absolute temperature 

T' (% - Dkj/CqrJ 

T-^ bulk temperature 

T b ^ T i " T b )k./(qr w ), E l uation ( 16 ) 

T^ absolute total temperature 

T t - T tV<<i r w ) 

T^ ^ bulk total temperature 

T t,b - T t,b) k i/(<l r w)> Equation (17) 

v r radial velocity 

v z axial velocity 



z axial distance from tube inlet 

z/r w 

0 p/°v 


% ^w/^z^bPb)^ shear stress parameter, 

Equation (23) 

£ second viscosity coefficient (£ = 0 for 

monatanic gas) 

(i viscosity 

P ‘ p/pi 

p density 

P * p/pi 

t wall shear stress 

T ' ^2,^1 : > 

Subscripts : 

b bulk, evaluated at bulk or mixed mean 

temperature 
i at tube inlet 

w at wall 


Superscripts : 

* on dimensionless quantity 
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Dimensionless axial distance, (z/D)/Re b 


Fig. 1. - Friction parameter based on local wall shear stress plotted against dimensionless distance from tube 
entrance. 



Dimensionless axial distance, (z/D)/Re b 


Fig. 2. - Friction parameter based on Ixal pressure gradient plotted against dimensionless distance from tube 
entrance. 


Nusselt number, Nu b » hD/l< b 





Dimensionless axial distance, z/D/Re b 


Fig. 3. - Friction parameter based on pressure drop plotted against dimensionless dis- 
tance from tube entrance. 



Reciprocal Graetz number, (z/D)/Re b Prj 

Fig. 4. - Local Nusselt number plotted against dimensionless distance from tube entrance. 
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Fig. 5. - Curves showing development of velocity profile. 
q‘ * 13. 5; Mi • 0. 06; Pr: - 0. 667. (Values of z/D are fpr 
Rej- 1500.) 



Fig. 6. - Effect of temperature ratio on Nusselt 
number for quasi-developed flow (Mj - 0). 
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Fig. 7. - Effect of temperature ratio on friction 
parameter for quasi-developed flows (M { * 0). 
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